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Abstract. We present an analysis of the rotational velocity of the primary of LB 3459 based on 107 new high-resolution 
and high-S/N ESO VLT UVES spectra. 105 of them cover a complete orbital period (0.26 d) of this binary system. We 
have determined an orbital period of P = 22 600.702±0.005 sec, a radial velocity amplitude of A\ = 39. 19±0.05 km/sec, 
and T = 2451917.152690*0.000005. From simulations of the He II A46S6A line profile (based on NLTE model atmosphere 
calculations), we derive v rot = 47±5 km/sec. 

We present an animation which shows the orbital movement of the binary system, its synthetic lightcurve, and compares the 
phase-dependent variation of the predicted with the observed He II A 4686A line profile. 

The radius of the cool component is almost the same size like Jupiter but its mass is about 70 times higher than Jupiter's mass. 
Thus, from its present mass (M2 = 0.066 M Q ), the secondary of LB 3459 lies formally within the brown-dwarf mass range 
(0.013 - 0.08 M ). It might be a former planet which has survived the previous common-envelope phase and even has gained 
mass. 
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1. Introduction 

LB 3459 is a close binary (P = 0.26 d) consisting of a sdOB pri- 
mary star and an unseen secondary with an extraordinary small 
mass. Based on the assumption that the primary component has 
a mass of Mi = 0.5 M Q , Hilditch et al. (1996, hereafter HHH) 
found that the cool secondary with M2 = 0.086 M is in ex- 
cellent agreement with lowest mass ZAMS models of Dorman 
etal. (1989). 

In a recent spectral analysis (Rauch 2000) of the primary, 
based on high-resolution CASPEC and IUE spectra, M\ = 
0.330 M Q has been derived from comparison of T e ff and logg 
to evolutionary tracks of post-RGB stars (Driebe etal. 1998). 
M2 = 0.066 M has subsequently been calculated from the sys- 
tem's mass function. Thus, the secondary is formally a brown 
dwarf(0.013-0.08M o ). 

HHH found that LB 3459 lies in that region of the 
Hertzsprung-Russell diagram appropriate for binaries under- 
going late case B (Roche lobe filled just before central he- 
lium ignition) mass transfer (Iben & Livio 1993) followed by a 
common-envelope (CE) phase. 




Send offprint requests to: T. Rauch 

* Based on observations made with ESO Telescopes at the Paranal 
Observatories under programme ID 66.D-1800. 
Correspondence to: Thomas.Rauch@sternwarte.uni-erlangen.de 



Fig. 1. The radial-velocity curve of LB 3459 (thick line). The phase 
is given with respect to the maximum positive velocity. The thin and 
dashed curves indicate the maximum velocity coverage during expo- 
sures of 180 and 3 600 sec, respectively, starting at a given phase. The 
phase coverage of one of our UVES spectra and of one CES spectrum 
used by Rauch (2000) are indicated. 

An alternative CE scenario has been studied by Livio & 
Soker (1984). They found that in star-planet systems, when 
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Fig. 2. Comparison of the co-added UVES spectrum to a CASPEC spectrum (thin line, observed by U. Heber 1985). Identified lines are marked. 



the star evolves into a red giant, the planet may survive the 
CE phase when its intial mass is above M cl -i t w 0.02 M Q . Then 
the planet can accrete from either the stellar wind or directly 
from the giant's envelope. Thus, it appears possible that the 
secondary of LB 3459 had formerly been a planet which has 
transformed into a low-mass stellar component. 

However, the mass-radius relation of LB 3459 derived 
from the mass function and the light-curve analysis (Kilkenny 
et al. 198 1) does not intersect with the empirical result of Rauch 
(2000). The reason for this is unclear. Two possible reasons 
are that the theoretical models by Driebe etal. (1998) are not 
applicable to the case of LB 3459 and the uncertainty for the 
photospheric parameters in Rauch (2000). 



Since the analysis of Rauch (2000) was hampered by the 
relatively long exposure times (lh and 2 - 3h, respectively) 
and hence, a relatively large orbital velocity coverage (the ob- 
served line profiles were broadened by the star's rotation as 
well as by a smearing due to orbital motion within the obser- 
vations), it has been speculated that log g from Rauch (2000) is 
somewhat too low. In order to make progress and to minimize 
the effects of orbital motion, spectra with short exposure times 
(180 sec) have been taken with UVES attached to the ESO 
VLT (Kueyen) which cover a complete orbital period (Sect. 2). 
From these spectra, the radial-velocity curve of LB 3459 is de- 
termined (Sect. 3). In Sect. 4 we will have a close look onto 
the spectral changes during the transit of the cool component. 
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Fig. 3. Radial- velocity curve of LB 3459 measured from the individual lines. He II A 4686A, Mg II A 4482A, and Si IV AA 3760, 4089A have the 
sharpest observed line profiles (shown here) while the Balmer lines are much broader and the S/N ratio is significantly decreasing towards the 
higher series members (Fig. 2). Note the velocity jumps close to T which are the result of the transit of the cool companion (Sect. 4)! 



Table 1. Orbital period of LB 3459 determined from different datasets. 
The typical error in the period is 0.05 sec. 



No. of measurements 


PI sec 


data source 


105 


22643.921 


Jan 8 


107 


22601.138 


Jan 6 + Jan 8 


131 


22600.702 


Jan 6 + Jan 8 + hhh 


24 


22590.326 


hhh 



From those spectra with a minimum effect of orbital smearing 
(v ra d at maximum, Fig. 1), we can then derive v rot (Sect. 5). 



2. Observations 

Two UVES spectra have been taken on Jan 6, 2001, with an 
exposure time of 180 sec each. Then the run was stopped al- 
though the night had been photometric. It was restarted on Jan 
8, 2001, and 105 UVES spectra were taken in a partly cloudy 
night. The improvement by the short exposure time is shown 
in Fig. 1. We used a slit width of l'/O projected at the sky and 
arrive at a resolving power R « 48 000. The spectra of Jan 8 
cover a complete orbital period of LB 3459. 

The spectra were subject to the standard reduction provided 
by ESO. The S/N ratio of the single spectra is about 20 - 30 
while from the results of the UVES exposure-time calculator 
we had expected S/N = 80 - 1 10 at a resolution of 35 000. In 
Fig. 2 we show our co-added spectrum: All 105 spectra (Jan 8) 
were shifted to the rest wavelength, then co-added and subse- 
quently binned to 0.045 A (3 wavelength points). Although the 
S/N of the co-added UVES spectrum is about 200, no chemi- 
cal element could be identified which had not been detected in 
other spectra before (Rauch 2000). 



3. Radial velocity 

In the first step, the phase-dependent radial velocity is de- 
termined by fitting Lorentzians to Hfi - H 10, He n /J4686A, 
Mg n A 4482A, and Si iv /U 3760, 4089 A in all 107 spectra. 
Other lines have been checked but turned out to be useless: 
C ii ^3921 A (too weak), Si iv X13921,4116A (blends with 
H e and C iv, respectively), and O in AA 3755, 3758, 3740A (too 
weak). During this procedure the smearing due to orbital mo- 
tion during the exposure was neglected because the exposure 
time (Sect. 1) is very short compared to the orbital period. 
Then, the derived velocity curves were fitted by sine curves 
in order to determine their periods, amplitudes, and Tq. 

The orbital period has been determined from the complete 
data set because the 107 spectra provide a time line of 9 periods 
which increases the accuracy. For the amplitude and r deter- 
minations, we use only the 105 spectra of Jan 8 because the 
scatter in the two spectra of Jan 6 increase the scatter in the re- 
sults artificially. The times (including heliocentric correction) 
are the middles of the exposures. The results are summarized 
in Fig. 3. 

Since the observed He n A 4686A line is the sharpest line 
in the optical spectrum (Rauch 2000), we achieve obviously 
(Fig. 3) the best fit to a theoretical sine curve with the veloc- 
ity curve measured from this line (the same line had been used 
by hhh). Mg ii /14482A is also sharp but much weaker and the 
Balmer lines are much broader and the S/N is lower. The mea- 
sured amplitude of the velocity curve appears to be larger for 
the higher Balmer lines. This is clearly an effect of the lower 
S/N. Thus, we decided to use only He n /14686A for the pe- 
riod determination. Moreover, hhh provide 24 additional mea- 
surements of the radial-velocity curve of LB 3459 based on 
He ii A 4686A which extends the time base up to 9692 peri- 
ods. In Table 1, we summarize the results based on different 
datasets. 
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Fig. 4. Comparison of radial-velocity measurements (Top: hhh data, 
bottom: our data) of LB 3459 with a sine curve calculated from our 
solution (P = 22 600.702 sec, A = 39.19km/sec, T = 51917.15269). 
Note that the measurements cover the time from Jan 29, 1994, to Jan 
8,2001. 

From the result of the dataset of Hen A 4686 A which 
includes additionally the hhh data (Table 1), we adopt 
P = 22 600.702 ±0.005 sec and Ai = 39.19*0.05 km/sec 
(Fig. 4). The period is only 0.02% longer than given by 
Kilkenny etal. (1991, 22 597.0 sec from photometric data 
of 27 eclipses). A higher radial-velocity amplitude of 
A\ = 40.8 ± 0.7km/sec has been measured by hhh from their 
observations of He n A 4686A. 

4. The transit of the cool companion 

The transit of the cool companion is clearly visible in the 
radial-velocity curves (Fig. 3). An inspection of the observed 
He ii A 4686A line profiles during this transit shows a decrease 
of the central depression close to To (Fig. 5). 

In order to simulate the phase-dependent variation of the 
line profile of He n A 4686 A, we have set up a surface grid (with 
25 azimuth and 25 altitude points) on the primary, and added 
up the spectra from these grid points with correct projected v ra d 
and consideration of limb darkening, calculated from the final 
model atmosphere for LB 3459 (Fig. 6). 

During the transit of the cool companion, the hidden parts 
of the primary's surface do not contribute to the resulting spec- 
trum (Fig. 8). It can be seen that the line profile during the tran- 
sit is first redshifted, then blueshifted, as can be expected. This 
can be seen as well in the radial-velocity curves (Fig. 3). At 
the begin and the end of the eclipse, the line profile appears 
deeper (Rossiter effect, Rossiter 1924). In addition, at the cen- 
tral eclipse, the line profile is shallower compared to phases 
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Fig. 5. The observed line profile of He II A 4686A in spectra taken dur- 
ing the transit of the secondary (times are given in HJD - 2451917). 
The spectra are binned to a resolution of 0. 1 A. The central depression 
is obviously smaller at T = 0.15196 (To = 0.15269). The dashed line 
indicates the rest wavelength, the vertical curve is the wavelength shift 
calculated from the radial-velocity curve (Sect. 3), the thin horizontal 
lines indicate the central depression outside the transit phase. 
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Fig. 6. Limb darkening in LB 3459 at four wavelengths as derived 
from our model atmosphere. Note that a standard darkening law (grey 
approximation) = 1+ ^° si:> (e.g. Unsold 1968, here calculated 

with /?o = 0.2, dashed line) is only a poor approximation for values 
which are calculated from a more realistic line-blanketed NLTE model 
atmosphere. 



shortly before and after. This is an effect of limb darkening, and 
it is qualitatively reproduced by our model atmosphere. From 
the comparison of the He n A 4686 A line profile at selected or- 
bital phases (tp = 355.7°, 359.3°, 6.5°), rxjr^ appears to be at 
its upper limit (Fig. 3). It is worthwhile to note, that a phase- 
resolved spectroscopic study of line-profile variation gives the 
radius ratio r\ / r2 independently from light-curve analysis — if 
the quality of the spectra and the phase resolution are sufficient. 

We have produced a short animation (available at the CDS) 
which shows all our observed spectra compared to the simula- 
tion as well as a view onto the orbital movement of the binary. 
In Fig. 7 we show a snapshot of this animation. 



T. Rauch and K. Werner: The rotational velocity of the sdOB primary of the eclipsing binary system LB 3459 (AA Dor) 



5 



. 342 344 346 348 350 352 354 355.7 358 359.3 2 4 6.5 8 10 12 14 16 18. 




_l I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 -4 -2 2 

wavelength / A 



Fig. 8. Top: Predicted phase-dependent variation (tp = 342°- 18°) of the theoretical He II .44686A line profile during the transit of the secondary. 
The dashed horizontal line indicates the central depression outside the transit phase. Bottom: Comparison of theoretical He II A 4686A line 
profiles at three different radius ratios, r\lr% = 2.50, 1.75, 1.13 with r\ = 0.168 ± 0.03 R G which are corresponding to the error limits of r\ and 
(Rauch 2000), with the observation (second panel, data binned to 0.1 A). tp e gives the angular limits where the primary is eclipsed. 




Fig. 7. Phase-dependent spectral variation of the He n A 4686A line 
profile. A snapshot taken from our animation. Left panels: View onto 
the binary system LB 3459 from the top and from the side (as seen 
from Earth, ( = 90°), at the bottom, the normalized lightcurve of the 
system is shown. Right: Our UVES spectra (bottom, smoothed by a 
Gaussian of 0.05 A FWHM) compared to the synthetic spectrum (see 
text for details). 

5. Rotational velocity of the sdOB primary 

The maximum change of v la d during a 180 sec-exposure of 
LB 3459 is Av rad = 1.958 km/sec (Fig. 1) and thus, the "line- 
broadening" effects of smearing due to orbital motion are 
small. Since the S/N of a single of our UVES spectra is about 
20 (in average), we decided to co-add all spectra (Fig. 2) in 
order to achieve the maximum S/N ratio. The alternative, to 
select only those spectra where Av la( j is below a certain limit, 
e.g. smaller than 0.5 km/sec, and to reduce the orbital smear- 



ing would reduce the percentage of spectra to co-add to 15 % 
— resulting in a much lower S/N. 

Since the resolution of the spectra is better than 0. 1 A (at 
the position of He 11 A 4686A), we checked whether it is neces- 
sary to consider the fine-structure splitting of He 11 A 4686A. In 
Fig. 9 we show that this results in a wider line profile and can- 
not be neglected, otherwise the helium abundance is strongly 
overestimated. 

We calculated a grid of H+He composed NLTE 
model atmospheres in the relevant parameter range, i.e. 
r eff = 40-45kK (in steps of 1 kK), logg = 4.95 - 5.55 (0.05, 
cgs), n H e/«H = 0.0001 -0.0020 (0.0001, number ratios), and 
calculated theoretical line profiles of He n A 4686A. We per- 
formed a x 1 test (wavelength, flux level, T e g, logg, «He/«H, 
and v ro t) in order to determine the rotational velocity. The low- 
est reduced^ 2 value (0.44) is found at r e ff = 44kK, log g = 5.4, 
«He/iH = 0.0007, and v rot = 43 km/sec. However, Rauch 
(2000) has determined T e ff = 42 * 1 kK from ionization equilibria 
(He i/He n, C m/C iv, N m/N iv/ N v, O iv/O v) which are 
very sensitive indicators for T e tf, and log g- 5. 21 ±0.1 from the 
broad hydrogen Balmer lines. Since the error ranges are small, 
we have to keep these values fixed. Then, v rot = 47 km/sec and 
«He/«H = 0.0008 yield the best fit. In Fig. 10 the co-added spec- 
trum is compared with theoretical line profiles of He n A 4686A 
which are convolved with rotational profiles with different ro- 
tational velocities and subsequently with a Gaussian of 0. 1 A 
(FWHM) in order to account for the instrumental profile. From 
this comparison we have estimated an error range of ±5 km/ sec 
for v r ot in this analysis. 

The «He/«H ratio is the same as found by Rauch (2000) who 
used CES spectra with a resolution of 0. 1 A. Since the exposure 
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Fig. 9. Theoretical line profiles of He II A 4686A calculated from mod- 
els with r eff = 42kK, logg = 5.2, and n He /n„ = 0.0008 (Rauch 2000). 
The thick lines show the profile which is calculated with fine structure 
splitting (components are indicated). For the thin lines, He II A 4686A 
is calculated as one multiplet. Note that the deviation is larger than the 
spectral resolution (0. 1 A) of our observation. At the bottom, the pro- 
files are convolved with a rotational profile (v rot = 45.7 km/sec) and 
compared to a theoretical line profile (dashed, without fine structure 
splitting) calculated from model with a 50% higher helium abundance 
(n He /n H = 0.0012) 
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Fig. 10. Comparison of theoretical line profiles (r cff = 42kK, 
log g = 5.2, n H Jn H = 0.0008) of He II A 4686A to the co-added UVES 
spectrum. At v rot = 47 km/sec we find a good agreement. The error 
range can be estimated to ±5 km/sec. 
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Fig. 11. Dimensions of LB 3459. 



time of the CES spectra had been 3600 sec, the orbital smearing 
(Fig. 1) and the lower S/N of about 20 made the fine-structure 
splitting (Fig. 9) unimportant. 

6. Dimensions of LB 3459 

With the orbital period P and the radial velocity amplitude A\ 
we can calculate the circumference U\ and radius a\ of the 
primary's orbit. From the relation a\M\ - CI2M2 we can derive 
parameters for the secondary. The components' separation is 
a — a\ + a%. Kilkenny et al. (1979) have given formulae for the 
radii of the components: r\ = 0.138a and r2 = 0.079a derived 
from the light curve. The calculated values are summarized in 
Table 2. In Fig. 1 1 the dimensions of LB 3459 are illustrated. 

We calculate r\ = 0.1678 ± 0.03 R for the primary. This 
is about 1/V2 smaller than r\ - 0.236 R G derived by Rauch 
(2000) from r\ = sjGMJg. Although Mj is determined by 
comparison of T e ff and logg to theoretical models (Sect. 1), g 
itself is the main error source in this formula — a change of 
0.3 in log g (factor 2 in g) would change M\ only by a few per- 
cent. Since all values of Table 2 have been calculated straight- 



forward, the errors are mainly propagation errors. This may be 
a hint for a higher g than determined by Rauch (2000). 

However, the velocity value of v rot = 45.7km/sec for 
bound rotation following Rauch (2000), i.e. calculated with 
r\ = 0.236 R Q and P — 22 600.702 sec would be smaller by 
1 / V2 as well. Then, the present rotational velocity is much 
higher than this value. For this paper, we decided to adopt 
r\ = 0.168 R and ri = 0.096 R©. However, the origin of this 
discrepancy needs further investigation. 

7. Results and discussion 

We have determined the orbital period of LB 3459 from the 
newly measured radial velocity curve in good agreement with 
previous values of Kilkenny etal. (1991). Together with data 
from hhh we could strongly reduce the error limits (Table 2). 

We derived a rotational velocity of the primary of LB 3459 
of v 10t = 47±5 km/sec which is in its error limits in agreement 
with v rot = 34±10 km/sec found by Rauch (2000). For r x = 
0.236 R , this is equivalent to a bound rotation. If the radius 
is r\ - 0.168 R Q (Sect. 6), then the primary rotates about 40% 
faster than bound. 



T. Rauch and K. Werner: The rotational velocity of the sdOB primary of the eclipsing binary system LB 3459 (AA Dor) 



7 



Table 2. Orbital parameters of LB 3459. For r\ and r 2 , see text for 
details. 
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The radius of the cool component is = 6.682 • 10 9 cm 
which is almost the same size like Jupiter (r = 7.143 • 10 9 cm) 
but its mass is about 1/15 M Q which is much higher than 
Jupiter's mass (1/1047 M Q ). Thus, from its present mass, the 
secondary of LB 3459 lies formally within the brown-dwarf 
mass range (0.013 - 0.08 M Q ). The numerical simulations of 
star-planet systems presented by Livio & Soker (1984) show 
that a critical mass exists above which it is possible that a planet 
with M 2 > M crit gains mass during a CE phase. Due to a sim- 
plified treatment of accretion and evaporation (Iben & Livio 
1993) in the work of Livio & Soker (1984), improved calcu- 
lations of star-planet systems are highly desirable to verify the 
idea of Rauch (2000), that the secondary of LB 3459 is a former 
planet. 

From the comparison of Hen A 4686 A to the co- 
added UVES spectrum (Fig. 10), the helium abundance of 
«He/«H = 0.0008 ±0.0002 found by Rauch (2000) has been con- 
firmed by the use of spectra with higher S/N and higher resolu- 
tion. Moreover, the rotational velocity found here is very close 
to the value of v rot = 45 km/ sec which had been adopted by 
Rauch (2000) in his analysis. 

Although the quality of the co-added UVES spectrum is 
much higher than previously available spectra, no elements 
have newly been detected in it. 

The reason for the discrepancy between primary and sec- 
ondary masses derived from spectral analysis and those from 
the radial-velocity and the eclipse curves (Sect. 1) is still un- 
clear — although some problems in the spectral analysis due 
to uncertain rotational broadening and smearing due to orbital 
motion could be ruled out. There is a hint, that the surface 
gravity g may be slightly higher than determined by Rauch 
(2000). If true, this would solve the problem with the discrep- 
ancy. Unfortunately, significant effects of changes in g are de- 
tectable only in the outer line wings where the data reduction 
of the broad Balmer lines in the echelle spectra is not very ac- 
curate. 



To summarize, there is no reason to perform a new spectral 
analysis of LB 3459 based on the new optical UVES spectra. 
However, it is highly desirable to obtain high-resolution and 
high-S/N spectra in the UV/FUV in order to improve the de- 
termination of r e (f based on ionization equilibrium of metals, 
to improve the determination of log g by the evaluation of the 
higher members of the Lyman series, and to identify other el- 
ements beside C, N, O, Si, Mg, Fe, and Ni in order to prove 
likely diffusion effects found by Rauch (2000). 
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